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In the context of the AdS/CFT correspondence we discuss the gravity dual of a high energy 
collision in a strongly coupled M = 4 SYM gauge theory. We suggest a setting in which two 
colliding objects are made of non-dynamical heavy quarks and antiquarks, which allows to treat the 
process in classical string approximation. Collision "debris" consist of closed as well as open strings. 
If the latter have ends on two outgoing charges, they are being "stretched" along the collision 
axes. We discuss motion in AdS of some simple objects first - massless and massive particles - 
and then focus on open strings. We study the latter in a considerable detail, concluding that they 
rapidly become "rectangular" in proper time -spatial rapidity r — y coordinates with well separated 
fragmentation part and a near- free- falling rapidity-independent central part. Assuming that in the 
collisions of "walls" of charges multiple stretching strings are created, we also consider the motion 
of a 3d stretching membrane. We then argue that a complete solution can be approximated by two 
different vacuum solutions of Einstein eqns, with matter membrane separating them. We identify 
one of this solution with Janik-Peschanski stretching black hole solution, and show that all objects 
approach its (retreating) horizon in an universal manner. 



PACS numbers: 

I. INTRODUCTION 



The AdS/CFT correspondence Q is a duality of the 
conformal (CFT) J\f—4 supersymmetric Yang-Mills the- 
ory and string theory in 5d Anti-de-Sitter space (AdS 5). 
Multiple papers use this fascinating theoretical tool, in 
a regime in which the gauge theory is in a strong cou- 
pling regime while string part is in weak coupling - the 
classical SUGRA regime. The equilibrium finite temper- 
ature version of this correspondence, using a black-hole 
background, was suggested by Witten Q. Applications 
of this version of correspondence to properties of strongly 
coupled high-T phase of QCD are very actively pursued: 
we will briefly review those in the next subsection. 

The aims of this series of works are however quite dif- 
ferent: instead of focusing on equilibrium thermal mat- 
ter, we hope to develop a gravity dual framework to time- 
dependent process of high energy collisions. We will not 
assume equilibration or use macroscopic variables like 
temperature or hydrodynamic flows: we hope to be able 
to understand how they naturally appear for collisions 
of large systems. Instead we focus on motion of strings 
in AdS$ in this work, and, in the second one, on "holo- 
grams" which an observer will see in our world - the 
AdS$ boundary - as a function of time. 

Since this is the first paper of the series, we decided to 
start with rather extensive introduction, which describes 
similar works and summaries our current understanding 
of the subject. 



A. Strongly coupled Quark-Gluon Plasma 

It is well known that non-perturbative properties of the 
QCD vacuum phase - confinement and chiral symmetry 
breaking - are absent above some critical temperature, 
where matter is in the so called Quark-Gluon Plasma 
(QGP) phase. Although at high T one naturally expects 
the QGP to be in a weakly coupled regime, it has been 
conjectured recently [|| that at least at T = (1 — 2)T C - 
known as the RHIC domain - it is closer to a 'strongly 
coupled' regime (sQGP). 

This was a significant "paradigm shift" in the field, 
and various directions toward the understanding of sQGP 
constitute a mainstream of the field. Basically there are 
two competing options: one, based on electric-magnetic 
duality [4|, relates small viscosity and diffusion of sQGP 
to presence of magnetic monopoles and predicts that it 
will disappear at T away from critical region. Another 
- based on AdS / CFT - relates it to "quasiconformal be- 
havior" of QGP at T > 2T C . A comparison between 
experimental results from RHIC (T = (1 — 2)T C ) with 
those at LHC (higher T) will hopefully shed light on it 
in near future. 

Let us only mention some important developments 
related to the latter approach, AdS/CFT. In a static 
finite-T setting with AdS-black hole metric Q the study 
started with classic results on bulk thermodynamics [5| 
and transport coefficients [6[: those works provided in- 
triguing results. It was shown that while the Equation 
of State can be quite close to that of weakly coupled 
plasma, the transport properties can differ from them by 
orders of magnitude. It is enough to mention that while 
viscosity to entropy ratio is believed to be limited by the 
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AdS/CFT value from below 



recent hydrodynamical studies by three groups Q have 
concluded that the experimental data on the so called 
elliptic flow are better reproduced if this ratio is even 
smaller than that! (For a possible way out of this puzzle, 
also based on AdS/CFT, see e.g. 0.) 

Then attention focused on high energy jet quenching, 
with the result that a heavy quark pulls a string, with 
specific and calculable shape. The AdS/CFT result for 
the drag force [ll| and heavy quark diffusion [l(| , turned 
out to be correctly related by the Einstein relation. For 
a recent brief summary see e.g. [l2j: it is sufficient to 
mention here that all these results seem to be in much 
better agreement with what is seen phenomenologically 
in heavy-ion collisions at RHIC than their weak-coupling 
counterparts. 

Further development of the jet quenching problem was 
related to the question where does the lost energy go?. 
In a hydrodynamical context it was suggested that the 
so called "conical flow" [HI of matter should develop, 
induced by a heavy charge moving in a strongly coupled 
plasma. The "hologram" of the dragging string has been 
calculated by Princeton and Seattle groups [3, EH]: it 
described the conical flow picture in stunning detail. 

B. Gravity dual for heavy Ion collisions 

The results mentioned above are all equilibrium ones, 
obtained using static AdS-black hole metric. Although 
they should be applicable for a macroscopically large and 
slowly expanding fireball, one may proceed to more de- 
manding issues related with AdS/CFT in time-dependent 
out-of-equilibrium setting. Those will provide new in- 
sights into equilibration issues, explaining when and with 
what accuracy thermo- and hydrodynamics become ap- 
plicable. 

In AdS / CFT language going from cold vacuum to hot 
plasma means going from pure AdS (extremal black hole) 
to black-hole AdS via creation of trapped surface. There- 
fore the problem to be considered is a kind of gravita- 
tional collapse, occuring in gravity-dual as a result of 
high energy collision. 

The quest for black hole formation in collisions has a 
long history we would not attempt to review here. Let 
us just mention that it was discussed for "real" grav- 
ity at colliders, which may get possible provided it gets 
strong due to extra dimensions. Black hole production 
in AdS spaces were discussed both in cosmological brane 
world models, as well as in AdS/CFT framework from 
late 90s: we only mention few papers most related to our 
work. Black hole emerging from collisions were discussed 
in AdS$ background by Horowitz and Itzhaki[l7|. who 
considered departing black hole. Giddings and Katz[!1| 
have discussed holograms of the falling objects in AdS$ 



background, in cosmological setting (which has some dif- 
ferences with AdS/CFT one in boundary conditions). In 
fill a solution for black hole creation from collision of 
particles was obtained for a simpler case of AdSs back- 
ground. It was recently further stuided by Kajantie et 
al. PH. 

In the context of gravity dual to heavy ion collisions, 
the problem of black hole formation was discussed by 
Sin, Shuryak and Zahed [2(| (SSZ below). One specific 
solution they discussed in the paper was a "hologram" 
of a departing black hole, corresponding to a spherically 
symmetric (Big-Bang-like) solution with a decreasing T 
. SSZ also proposed two other idealized settings, with d- 
dimensional stretching, corresponding for d=l to a colli- 
sion of two infinite thin walls and subsequent "Bjorken" 
rapidity- independent expansion [21(, with 2d and 3d cor- 
responding to cylindrical and spherical relativistic col- 
lapsing walls. 

Janik and Peschanski 22] (below referred to as JP) 
have addressed the simplest wall-on- wall collision. In this 
case the time and longitudinal coordinate x\ are natu- 
rally substituted by the proper time and spatial rapidity 

r=^ 1 ,y= 1 -lo 9 (i^) (2) 

since the rapidity-independent solution depends on only 
r. Instead of solving Einstein equations with certain 
source, describing gravitationally collapsing "debries" of 
the collision, JP applied an "inverse logic", extrapolat- 
ing into the bulk the metric which yield expected hy- 
drodynamical solution at the boundary. JP found an 
asymptotic (large-time) solution for a "stretching AdS- 
BH" . As expected, it indeed possesses a horizon moving 
away from the AdS boundary, as Zhorizon ~ t 1 / 3 . A very 
important feature of the leading-order JP solution is en- 
tropy conservation: is that while their presumed horizon 
is stretching in one direction and contracting in others, 
to the leading order two effects compensate each other 
and keep the total horizon area constant. We will discuss 
a bit more this solution and use it in section HV Al 

Further discussion of the subleading (next power of 
inverse time) terms has been made by Sin and Nakamura 
[23[ who identified corrections to the JP solution with 
the viscosity effects. Terms of still higher order have 
been subsequently studied [24[ , but eventually Janik et al 
[25l ] concluded that the expansion series are inconsistent 
beyond the first few orders. Our view is that this is 
how it should be, and the arising near-horizon singularity 
indicate that presence of matter term (absent in JP) is 
inevitable. 

Unlike JP et al we will not use any "inverse logic" and 
will not be looking for the solutions corresponding to pre- 
determined hydro on the boundary. Instead we will focus 
on the formation stage, whether black hole is or is not 
formed, and will calculate the (time-dependent) stress 
tensor on the boundary, whether it is a hydro-type on 
not. 
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C. Hadron collisions in QCD, the Lund model and 
the "Color Glass" 



Rather early in development of QCD, when the notion 
of confinement and electric flux tubes - known also as 
the QCD strings - were invented in 1970's, B.Andersen 
and collaborators [26[ developed what gets to be known 
as the Lund model of hadronic collisions. Its main idea is 
that during short time of passage of one hadron through 
another, the strings can get reconnected, and therefore 
with certain probability some strings become connected 
to color charges in two different hadrons. Those strings 
get stretched longitudinally and then break up into parts, 
making secondary mesons and (with smaller probability) 
baryons. Many variants of string-based models were de- 
veloped, and some descendants -e.g. PYTHIA - remain 
widely "event generators" till today. 

If there are several string stretched, it is usually as- 
sumed that both their interaction and influence on break- 
ing is negligible. 

However if one either considers very high energy col- 
lisions, when a single hadron should be viewed as being 
made of many color charges (partons), or heavy ion col- 
lision, a different asymptotic picture has been proposed. 
McLerran and Venugopalan [28| argued that instead of 
multiple string the fields produced should be considered 
as classical gauge fields -known as Color Glass model - 
and their subsequent evolution be derived from solution 
of classical Yang-Mills equation [29| . They suggested this 
regime is true at very high parton density, when the ef- 
fective coupling is weak. Accepting the Color Glass pic- 
ture as a correct asymptotic for very high parton density 
and large saturation scale Q s —> oo, one still wanders 
what should happen in the case of intermediate scale 
Q s ~ .3- 1.5 GeV. 

Recent developments of the so called AdS/QCD pro- 
posed a view that this interval of scales in QCD constitute 
a "strong coupling window" . In particular, Brodsky and 
Teramond [27j have argued that the power scaling ob- 
served for large number of exclusive processes is not due 
to perturbative QCD (as suggested originally in 1970's) 
but to a strong coupling regime with near-constant cou- 
pling (quasi-conformal regime). Polchinski and Strassler 
[3(J | have shown that in spite of exponential string ampli- 
tudes one does get power laws scaling for exclusive pro- 
cesses, due to convolution (integration over the z vari- 
able) with the power tails of hadronic wave functions. 
One of us proposed a scenario [3l[ for AdS/QCD in which 
there are two domains, with weak and strong coupling. 
The gauge coupling rapidly rises at the "domain wall" as- 
sociated with instantons. Such approach looks now nat- 
ural in comparison to what happens in heavy ion/finite 
T QCD, where we do know that at comparable parton 
densities the system indeed is in a strong coupling regime. 




FIG. 1: Schematic view of the collision setting. The classical 
heavy charges move along directions x± and collide at the ori- 
gin. String snapping leads to longitudinally stretched strings 
(wide black line) which are also extended into the 5-th coor- 
dinate r toward the AdS center at r = 0. The heavy charges 
move on the plane r = oo 



D. The goals of this series of papers 

In short, it is to study self-consistently the collision 
process in AdS/CFT. For hadronic collisions we basi- 
cally follow QCD-string-inspired (Lund) picture of the 
collision. While QCD phenomenology focused on "string 
breaking" , in AdS/CFT setting we will have instead their 
"falling" (departure from the boundary) into the IR. 

In this paper we will study in detail motion of "de- 
bris" - massless and massive particles and open strings, 
and membranes - in AdS$. In the second paper we will 
calculate the corresponding "holograms" of these objects 
- the stress tensor of matter created on the boundary. 
Although "debris" fly away into the 5-th direction, the 
usual energy and momenta are conserved in our world, 
and those "holograms" describe a flow of matter out- 
ward from the collision point. As we mentioned already, 
this can be viewed as a strongly-coupled version of Color 
Glass, put in the realm of A/"=4 SYM theory. 

We hope in subsequent works to go beyond the lin- 
earized gravity and follow nonlinear effects leading to a 
gravitational collapse of debris and formation of trapped 
surfaces. This would be dual to information loss (entropy 
production) and appearance of equilibration. 



II. THE SETTING 

One important suggestion made by SSZ is that heavy 
ion collisions posses "some internal high momentum 
scale" , usually called Q saturation , related to high density 
of color charges in boosted heavy ions. In order to model 
it more simply, we now propose substitute energetic light 
quarks by heavy ones, with the mass Mq of heavy fun- 
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damental quarks Q introduced into AdS/CFT via D7 
brane [32I]. As soon as Mq is at the scale of Q saturation, 
it makes little dynamical difference: but in the AdS/CFT 
language treatment of heavy quarks is simpler, as they 
are sources of classical strings. (This simplifying feature 
has been put to heavy use in treatment of the heavy 
quark jet quenching 

We will further assume that heavy quarks have no dy- 
namics of their own, as they are moving along straight 
lines 

x± = X\ ± vt (3) 

with constant velocity v, both before and after the col- 
lisions, see FigfT] If so, there is no conventional gluonic 
radiation on the brane or gravitational radiation from 
them in the bulk, as there is no acceleration. 

The dynamical objects we will focus on are classical 
strings, ending at these heavy quarks and propagating 
in the bulk (for metrics changing from AdS to JP-like 
one). We will study which solutions exist as a function 
of collision rapidity and whether they are stable or not: 
we will conclude that at sufficiently large v > v c these 
strings basically go into free fall toward the AdS center. 

The next step is to consider not a single pair of charges 
(a single stretching string) , but many. One limit is a pair 
of colliding "walls of matter" , containing multiple heavy 
quarks. For simplicity, think of these two walls as CP 
mirror images of each other, made of colorless "dipoles" . 
"Snapping" of their string at the collision leads to mul- 
tiple strings, all of which being stretched longitudinally. 

We then argue that many such strings combined could 
be considered as a thin singular sheet of matter, referred 
to below as "membrane" . (Note an important distinction 
between a membrane and a "true brane" : since the for- 
mer has only energy-momentum but lacks the RR charges 
and consequent Coulomb repulsion, it cannot "levitate" 
like branes, and simply falls under gravity.) 

It has been shown by Israel [33[ how a gravitational col- 
lapse of a thin layer of matter can be described via two 
different discontinuous vacuum solutions of the Einstein 
equation without matter (T^ — 0). Self-consistency of 
the solution is then reached by fulfilling covariant junc- 
tion conditions, resulting in membrane equation of mo- 
tion. 

The issue of self-consistency will not be addressed in 
this work: we will discuss below falling of various objects 
- particles and open strings, as well as 3+1 membranes - 
ignoring for now the effect of their own weight on the met- 
ric. The proposed evolution of the system is explained 
schematically in Fi^5] Part (a) of it shows some snap- 
shots of this surface, at some early time and then at a 
later stage. The horizontal direction is the collision di- 
rection x\ while the one along the circles represent any 
of the two other transverse directions x%, X3 (on which no 
dependence is expected). The radial direction r in part 
(b) is the 5-th AdS radial direction, a distance from the 
AdS center. Since the "membrane" is being stretched in 
x\ (linearly in time), it has to retreat in r and become a 



thinner cylinder, just as a stretching soap film will do in 
a similar setting. 

At this point we would like to emphasize a close anal- 
ogy, as well as differences, with the jet quenching prob- 
lem. One studied first a single falling string governed by 
simple Nambu-Goto action and the overall metric. The 
complicated picture of matter flow is then recovered us- 
ing weak (linearized) gravity. One difference is that in 
a jet quenching problem the string is stationary (in the 
charge frame) while in our case it is not. Furthermore, 
we will discuss also multiple strings, which may form an- 
other singular object - the membrane. Also the metric 
in our problem is first considered to be just AdS, but 
eventually it will be non-trivially affected by the mem- 
brane's own weight. If so, one should no longer use the 
linearized gravity but solve Einstein equations in its full 
nonlinear form. 

Needless to say, this is a very difficult task, amenable 
to analytic treatment only if some drastic simplifications 
are made. A scenario outlined in Fi^2fa) would have met- 
ric dependent on 3 variables: time, longitudinal direction 
and the AdS radial one, t,x\,r. We thus propose a fur- 
ther simplification of the problem: changing variables to 
proper time and spatial rapidity |(5J) we would look for 
y-independent solutions, corresponding to purely cylin- 
drical part of the membrane in the middle of Fij^a), 
ignoring the curved "fragmentation" regions. With only 
two variables, r, r one has a problem of similar level of 
complexity as the one addressed by Israel [3^|, for a spher- 
ical gravitational collapse. 

Further clarification of the proposed scenario is shown 
in Fig(2)(b), displaying a trajectory of the membrane r(r). 
During the first stage of the process the "debris" of a 
collision in a bulk - the particles and open strings - are 
accelerated by the AdS gravity and fall into the 5-th di- 
mension till they reach the relativistic velocity v « 1 
(stage 2)). If there be only one object falling, its grav- 
ity being negligible compared to overall gravity of the N 
branes at the AdS center and they would simply continue 
their relativistic fall. However large number of them have 
enough mass to create a horizon which suddenly slows 
down the membrane (as a distance observer sees it[4(| ): 
at stage 3 the membrane is trailing the receding horizon 
(the dashed line). 

If we would discuss pure AdS/CFT theory this would 
be the end of the story: but in other more QCD-like set- 
ting one can have an additional potential which will stop 
membrane because of existence of a stationary "decon- 
finement" horizon. If so, the system reaches a "mixed 
phase" era with stationary horizon and fixed T, similar 
to static fireball discussed by Aharoni et al [3J] except 
that in our setting the longitudinal stretching continues. 

The trajectory of the collapsing matter sheet should 
be such as to provide a consistent solution to Einstein 
equations, combining the JP-like vacuum solution outside 
the falling sheet, with the "stretching AdS" inside it. 

The paper is structured as follows. In the next section 
we solve equation of motion for different objects falling 
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FIG. 2: (a=upper) Two snapshot of the membrane shapes, 
at different time moments. See text for explanation of the 
coordinates. (b=lower) Schematic view of the four periods in 
gravity dual solution in which falling objects are (1) acceler- 
ated into the 5-th dimension r till they reach a relativistic 
velocity v « 1, then (2) continue their relativistic fall till (3) 
breaking near the retreating horizon. 



in AdS. We start with massless and massive particles in 
subsection IIII Al 

The main part of this work is study of the open strings, 
being stretched between two departing charges. We de- 
rive analytically the so called scaling (factorisable) solu- 
tion in section IIII Bl Similar solutions have been used 
previously in connection to anomalous dimensions of 
"kinks" . New part is discussion of the limits for its exis- 
tence and stability. 

We then find more general non-factorisable solutions 
in section Till CI which can only be obtained numerically 
We find that in proper time -spatial rapidity coordinates 
r, y we use those basically becomes "rectangular" , with 
a nearly free- falling rapidity-independent part. We con- 
clude this section with results for falling membranes. The 
next section starts with an introduction to the issue of 
"stretching black holes" in section IIV Al and concludes 
with section HVBI in which we show that all objects con- 
sidered above are approaching the (retreating) horizon in 
a very universal fashion. We conclude with some discus- 
sion and outlook in section [Vj 

In the second paper of the series we will calculate back 
reaction of gravity, by solving linearized Einstein equa- 
tions and obtaining stress tensor on the boundary ( "holo- 
grams") for some of these falling objects. 



III. OBJECTS FALLING IN AdS 5 

The collisions creates a lot of "debris" in form of vari- 
ous excitations. Since we would like to follow the collision 
in the bulk, we naturally have to think of them in terms 



of string theory. Thus there are the following types of ob- 
jects: (i) massless and massive particles; (ii) open strings, 
with ends at the receding walls; (iii) membrane. The 
"open string" category is naturally split into "mesons" 
with both ends on the same wall, and "stretched strings" , 
with both ends attached to different walls and moving in 
the opposite direction. We will consider a set of multiple 
strings copied many times in transverse dimensions X2 , X3 
as a 3-d membrane. The validity of this approximation 
will be explained later. 



A. Falling particles 

As is usually done in this kind of problems, the AdS 
radius is inverted, so that a coordinate z = 1/r is used 
instead of r. The AdS boundary is thus at z — and 
"falling" objects move away from it toward infinity. The 
AdS§ x S 5 metric in such coordinates is 



ds 2 



-^-{dx 

z A 



dt 2 



dz 2 ) 



r 2 <m\ 



(4) 



where the last term, related to angles of S 5 is of no impor- 
tance in this work. We choose to work in t, y coordinates 
mentioned above ©. The metric is translated into the 
following form: 



ds 2 



R "-(-dr 2 +T 2 dy 2 + dz 2 ) 



(5) 



where we ignore the transverse coordinates and the S$ 
part. 

One feature of AdS§ metric is its boost invariance, the 
importance of which will be seen later. Let us assume 
particles move with constant spatial rapidity y, so the 
trajectory can be described by z{t). Massless particles 
move along the geodesies with zero interval ds 2 = which 
in the metric ([5]) simply means z = r. 

Massive falling objects were already discussed in [37| . 
but here we present it in a different form, more closely re- 
sembling much more nontrivial ones in the next sections. 
Using the coordinate time r one simply write down the 
interval as an action for a particle moving in the 5-th 
direction of 



S 



dr 



z(r) 



(6) 



where the non-trivial trace of the AdS metric is z in the 
denominator. This leads to well known EOM 



z(r) 



z(r) 



(7) 



Nonrelativistically, one can neglect i(r) and think thus 
about a motion in a logarithmic potential well[4l|. Ul- 
trarelativistically, one finds instead that as z(r) — > f the 
acceleration goes to zero, as needed. Thus, in the stan- 
dard coordinates, very little seems to happen after the 
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particle reaches ultrarelativistic regime: it runs forever 
toward z — > oo with speed of light. But this is a (well 
known) illusion due to relativistic time slowing: in its 
own proper time, the particle continue to accelerate and 
reaches the AdS center in finite proper time. 
This EOM is easily integrated yielding 



z(t) = \Jt 2 + VqZqT + Zq 



(8) 



B. Falling open strings: the scaling solution 

After this little warm-up, let us consider motion of the 
open strings. Its action is given by Nambu-Goto, and if 
one ignores two transverse coordinates xi , X3 and uses as 
two internal coordinates the t, x (time and longitudinal 
coordinate) the string is described by by one function of 
two variables z(x,t). The corresponding string action is 
then 



5' 



R 2 
2vra' 



dt 



dx 



K dx } K dt> 



(9) 



Note that only one term, the time derivative, is differ- 
ent from long-used static action used in [35] for static 
calculation of the inter-charge potential. The boundary 
conditions would be z = at two rays x = ±vt, the 
world lines of the heavy quarks. (The boost invariance of 
the AdS$ metric allows us to work in a frame where the 
open string endpoints move with opposite velocities) 

Translating into the r, y language, the boundary con- 
ditions are now determined at fixed y = ±Y where 
v = tanhY and Y is the rapidity of the heavy quarks 
(colliding walls), by doing so, we transfer time depen- 
dence from the boundary conditions into the equations 
themselves. The corresponding action is now 



S = 



R 2 



rdrdy 



V 



dz 



(10) 



Before solving the corresponding equation in full, we 
will first discuss "scaling" solutions in the separable form 



(ii) 



suggested by conformal properties of the theory. Such 
solutions were known in literature [36[ , in Euclidean con- 
text, they were used for AdS/CFT calculation of the 
anomalous dimensions of "kinks" on the Wilson lines (of 
which our produced pair of charges is one). 
The scaling ansatz leads to a simple action 



S = 



_R_ 

2na' 



drdy 



V/' 2 + / 4 - / 2 



(12) 



Using the fact that y does not appear in the action, there 
is a conserved "energy" 



V 



VF + v 



= E 



(13) 
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FIG. 3: Rapidity of the collision Y = arctanh{v) vs f§ . The 
maximum gives a critical rapidity Y c . For Y < Y c , two f§ are 
possible, corresponding two string configurations. For Y = Y c 
,only one /q is possible. The region Y > Y c cannot be reached 



with the "potential" V = f 4 — f 2 , and thus the derivative 
of the function / can be readily obtained 



/' = 



VV(V-E 2 ) 
E 



(14) 



Note that the function / decreases from infinity on the 
boundaries to its lowest value at the middle of the string 
which we will call /o, so / > /o- At / = /o the derivative 
vanishes, so (|T4"|) provides also a simple equation /q — 
/ 2 - E 2 = relating E to / . 

Integration of (fl4|) gives the following solution 



V = fo\ 
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15) 



where F and 77 are elliptic integral of the first and 
the third kind. f§ depends on collision rapidity Y = 
arctanh(v) via the boundary condition at f(Y) — 00, as 
shown in Fig. [3] 

The existence of a maximum means that there are no 
scaling solutions when the rapidity Y is larger than some 
critical value, while if the quarks move on the boundary 
slower that the critical rapidity, there are two solutions. 

In order to characterize the solutions, it is useful to 
introduce "effective potential" for two separating quarks 
for each scaling solution, defined as instantaneous energy 
U = AS /At, where AS is action given by the area of 
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the string world sheet, At is the time interval. U needs 
to be regulated, which is obtained by subtracting the 
Wilson loop corresponding to two non-interacting moving 
quarks. In other words, we calculate the subtracted area: 




df (16) 



The second term corresponds to /' = oo, precisely the 
straight string going in z direction, which is AdS solution 
for a moving quark. Note that we have switched to t,y 
coordinates, which does not change the form of the string 
action (fT2"|) . With this prescription, we calculated U for 
solutions in both branches, which are compared in Fig. [4] 
The solution with the lower potential has a chance to be 
the stable one, while the higher potential one (with large 
fo, or longer string) must be metastable. 

Let us now comment on the small v limit of the scaling 
solution. At large separation (realized at late time) the 
quarks can be considered as quasi-static. At small v, or 
large f§, the effective potential can be simplified to the 
following form 




i? 2 / 

= -tj— t -0.5991 y/f - 0.1780 - 
2ira' \ jo, 

and relate more simply the velocity and fo 

0.5991 0.03115 



1 \/t 



fo 



fo 



(17) 



(18) 



Combining (jTTJ) and (fl"8"|) , we obtain the effective po- 
tential for small velocity and large separation to be 



V = 0.2285 



'1 + 0.6830 v 2 ) y / g 2 N 



(19) 



The coefficient in front (the leading term at v —* 0) coin- 
cides with the well known coefficient of static Maldaccna 
potential. 

The second term is thus the velocity-dependent "Am- 
pere's law" 0(v 2 ) correction to it. We are not aware of 
any other previous calculation of this term, except for 
the paper by Zahed and one of us 38] in which, based on 
resummation of ladder diagrams via Bethe-Salpeter eqn, 
the result was that the velocity dependence is 



U(v)/U(v = 0) = VI - viv 2 fa 1 + .bv 2 + 



(20) 



It is close but not the same[42j. 

Both branches of the scaling solution was also con- 
firmed by solving the equation numerically, starting from 
the middle point and scanning all values of fo- 
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FIG. 4: the potential V as a function of v for different 
branches of solution, circles for large- /q bran ch, crosses for 
small- fo branch V is plotted in unit of a/ 'g 2 N / L The poten- 
tial from the large /q branch is lower than that from small /□ 
branch 



The applicability of the scaling solution for a particu- 
lar Y depends of course not only on availability of a so- 
lution, but also on its stability i.e. how does the scaling 
solution evolve with time(r), given some perturbation at 
initial time. Denoting scaling solution g s (y) 
perturbation as 



-77-T and 

f(v) 



z{r,y) = Tg(r,y) 



](T,y)=9s{y)+6g(T,y) (21) 



we want to know whether the perturbation will grow or 
decay with time. The EOM for g(r,y) 



ay or ay dray 



■9g, 2 d 2 9 , o d 2 g 2 dg 



d 2 g dg^ 



T 9{^r) Tj-o + 2r— g — + r ^g(^-) - 2— - g 



y dy dr 2 dy 2i> dr 

d 2 9 2 , 7 d 9 o 3 d 9 
■-5-59 +7rs— -3rg — 
ay z or or 



dy 2 x dr J 
' dr' 



3tV(^) 2 + - 2 / 5 



dy 
&\ 
'!h 2 



oy or OTOy 



can be used by plugging (|2"Tj) in (f2"2"|) , and keeping only 
term linear in 5g{r 1 y) (consider only sufficient small per- 
turbation), we obtain the following linearized EOM for 
the perturbation: 



or oy 

-- 



D 



d 2 



drdy 



E 



dr 2 



F 



dy 2 



fig(T,y) 

(23) 
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with 



A = g"gl + 6. 9s - 



B = r(2g s g' 2 

C = -4g' s 

D = -r(2g 2 g' s ) 
E = r 2 (g s g' s 2 + g s ) 
F = gl - g s 



g s 

7g s - 3 5s 3 ) 



(24) 



define f = Inr as our time, the EOM simplifies to: 



= 
with 



4 

OT 



4 

ay 



b 32 



dfdy 



E W 2 



dy 2 



Sg(r, y) 
(25) 



A = A,B = B-E,C = C,D = D,E = E,F = F (26) 

(To make it easier to get all these functions one can ap- 
proximate scaling solution g s {y) with some parameteri- 
zations: we found that (^jj) 3 + (y )" = 1 fits all the 
scaling solution very well.) 

We need to seek eigenfunction Sg(r,y) = e Xr ip(y) sat- 



isfying (j2"5|) and boundary condition ip(y — ±Y) = In 
general, out of many eigenvalues A we should be inter- 
ested in those with positive real part, which will allow us 
to conclude when the solution is unstable. 

The eigenfunction results in the following EOM: 



C 



d d 2 
Gin — I - Ci-^-R 
ay oy z 



<P(y) = o 



(27) 



with 



Co = X 2 gs(g? + 1) + X(g s gf + 6g s + 2g 2 s g" s - 3g 3 s ) 

+3g'a9a + Q 9s - Ws - g"s 

C\ = -2g' s (\g 2 + 2) 
C 2 =g s (g 2 -1) 

Due to the symmetry y «-> — y of the problem, we can 
solve it in the positive-?/ region, with boundary condition 
ip(Y) — 0,t/)'(0) = 0. To solve this Schrodinger-like eqn, 
we use the iterative method. Starting on one boundary 
with tjj'(0) = 0,-0(0) = 1, the second condition only af- 
fects the normalization of ip(y). With some initial value 
of A, we can obtain the ip(Y) from the EOM. then we 
variate the value such that tp(Y) converge to 0. The re- 
sulting A gives the eigenvalue. Without much difficulty, 
we found the following set of eigenvalue for different Y, 
shown in TableHJ We also plot the eigenvalue A in the 
complex plane Fig[5] The evolution trend of this set of 
eigenvalues suggests that the transition from stable to 
unstable occurs at Y m inside .22-. 27 interval, which is 
way below the critical value Y ~ .5 above which there 



TABLE I: one set of eigenvalue for different rapidity 



A(i(n 2 ) 


4.2+94.8i 


3.3+I26.7i 


2.8+157.5i 


2.0+I88.5i 


Y 


0.48 


0.45 


0.42 


0.40 


\(w-' z ) 


1. 2+222. fi 


0.78+265.7i 


0.38+299.5i 


O.I2+346.4i 


Y 


0.37 


0.33 


0.30 


0.27 


AtKT 2 ) 


-0.27+404.2i 


-0.63+492.9i 


-0.80+569.8i 




Y 


0.24 


0.21 


0.18 





500 



400 



Im(lambda) 



300 



200 



100 




12 3 4 

Re(lambda) 

FIG. 5: The evolution of eigenvalue A from Y=0.48 to 0.18 in 
the complex plane 



were no scaling solutions at all. This shows that we essen- 
tially lose the scaling solution to instability for Y > Y m : 
we were not able to tighten this limits any further. 

In summary, the scaling solution exist only for suffi- 
ciently small rapidities Y < Y c ~ 0.5. Furthermore, we 
were able to verify that it is classically unstable already 
for Y > Y m w 1/4. Therefore solutions other than the 
scaling one is need for large rapidity, which is more im- 
portant for our purpose. 



C. Falling strings: the non-scaling solutions 

In this section we study generic solutions outside the 
scaling ansatz. But before we do so, let us explain quali- 
tatively why such solution must fail as the rapidity of the 
collision grows. The scaling solution, in which t and y 
dependences factorize, means that one tries to enforce a 
particular stable profile to a string. But as the rapidity 
gap 2Y between the walls grows, we so-to-say try to build 
wider and wider "suspension bridge" out of the string: it 
is going to break under its weight at some point. 

We again use z(r,y) = Tg(r,y) and EOM The 
boundary condition is g(r, y — ±V) = 0. Due to the sym- 
metry of the problem, it is sufficient to solve the dynam- 
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FIG. 6: The dynamics of the string(half ) g(r, y) with y — 0.6. 
The profiles from the innermost to the outermost correspond 
to t = l(solid red),r = 2(dotted blue),r = 4(dashed green), 
r = 8(dot-dashed black). 



ics of half of the string, with initial condition g(r, Y) = 
and|(r,0) = 0. 

However there are two potential problems in ([22]) . 
(i)the y derivative diverges on the boundary, (ii)the PDE 
is highly nonlinear and will show self-focusing of energy 
at certain "corners", as we will see. These make it dif- 
ficult to obtain a well-behaved numerical solution [43| . 
and to improve the performance of Maple PDE solver 
we used function h(r, y) = g(r, y) n as dynamical vari- 
able, with properly chosen integer power n so that the y 
derivative is finite on the boundary. 

Fig[6] shows the dynamics of the string with Y = 0.6. 
We start from the initial condition (^ygf^) 3 
and gfr(l, y) — 0. We choose the initial time r 
avoid the singularity at r = 0. n = 6 is used in solving 
the PDE. As time grows, the string profile approach a 
rectangular shape with sharper and sharper turn at the 
"corners" . Based on the numerical solution, we infer that 
in the r, y coordinates, any point of the string other than 
the boundary will ultimately become free falling when 
time is sufficient large. This can be supported by the 
following qualitative argument. Any tiny piece of string 
experiences the AdS effective gravity and the drag from 
its neighbors. Since in the non-scaling solution, the whole 
string keeps falling, it is natural to expect any point of 
the string approach the speed of light asymptotically, end 
up with a rectangular profile. Therefore, we conclude the 
edge of the profile is not important asymptotically. It can 
be well approximated by a flat profile in y, which will be 
studied in the next section. 



(#) 



3 = 1 
1 to 



D. Falling strings and membrane in AdSa 

The falling string can be considered as a solution at the 
center of the generic case considered above in the large 
rapidity limit of the ends Y — > oo. which makes z y- 
independent. Ignoring all derivatives over y in the EOM 
above one gets an ODE problem with the following eqn: 



-2r 



zz — z 3 z + T ZZ + 2 z 2 t = 



(28) 



which is similar but not identical to that of a falling mas- 
sive object ([7]): the difference comes from dimensionality 
of the object: 1/z 2 in the action (instead of 1/z), be- 
cause the string action is a 2-dimensional integral. It is 
now explicitly depending on r: there is no integral of 
motion but one can straightforwardly solve the EOM for 
different initial conditions numerically. We found at large 
r, g tends to 1. Therefore we show in this extreme case 
that the asymptotic solution is again z ~ t 

Summarizing the falling of all string objects, they have 
a universal asymptotic behavior z ~ r. Therefore we may 
model the falling particles/open strings by a membrane, 
which is made of multiple strings and is flat in 23 and 
y coordinates 

The coefficient in its DBI action, the membrane ten- 
sion, is now proportional to the density of charges in 
the colliding walls, and thus can be very large. This fact 
would mean that the membrane should eventually be con- 
sidered heavy enough, so that its weight would affect the 
metric itself. Since in this work we would not attempt 
to solve this problem yet, we treat the membrane as a 
test body falling in external AdS metric. In this case the 
value of its tension does not matter, and the action is 
very similar to Nambu-Goto string action except of the 
different power of z (now 1/z 4 ) 



Tdrdydxidx; 



(29) 



We parametrize the membrane with r,y,X2,Xa, and z- 
coordinate is a function of r only, z = z(t). The EOM is 
readily obtained, it is similar to the y-independent string 
case (coefficients 2 change to 4 in two terms): 

zz - z 3 z + Atz 2 + t'zz - At = (30) 

Its asymptotic solution is again z ~ r. 



IV. NEAR-HORIZON "BRAKING" 

A. Stretching black holes 

The JP solution we will now discuss addresses the first 
case, d=l. The main feature of the JP solution is that 
these two variables enter the metric via one specific com- 
bination 



v 



(31) 
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which simplifies Einstein's eqns and leads to a solution. 
JP have found that only for one particular power 7 = 
1 /3 there is no singularity at the horizon in one of the 
invariants - the square of the 4-indcx Riemann curvature, 
and argued that thus this solution should be preferred on 
this ground. 

However it is not clear what the physical meaning 
and significance of this singularity may be, in general. 
Furthermore, in the "membrane scenario" proposed in 
this work the JP-like metric only extends from the AdS 
boundary till the falling membrane, while the would-be 
singularity is in the second domain, where this solution is 
not supposed to be used at all. It is, so to say, a "mirage 
behind the mirror", singular or not does not matter. 

There is another reason why this particular power 
should be selected: only in 7 = 1/3 case such that the 
total area of the horizon (3d object normal to time and 
z) is time independent: the factor r (from stretching y\) 
is canceled by the factor 1/z 3 from contracting z. Thus, 
this stretching solution is area-preserving, and thus po- 
tentially dual to the entropy-conserving adiabatically ex- 
panding fireball. 

The specific form of the JP metric is 



had canceled out. Also note that there is a significant 
singularity at the horizon (v — 3 in these units) in this 
stress tensor, which is again irrelevant because this metric 
is not supposed to be used there. 



B. Objects approaching the horizon 

Before we discuss the JP metric, let us remind the 
reader how this approach works in the usual black holes 
with the Schwartzschild metric: it will be needed to em- 
phasize the difference between them. 

Massless particle falling radially in the Schwartzschild 
metric satisfies the ds 2 = eqn, which is 



y dt' 



(1 



r 



leading to exponentially fast "freezeout" , 

(r - r h ) ~ exp(-t/r h ) 

The same is also true for other objects, of course. 
We use the following rescaled coordinates: 



(37) 



(38) 



ds 2 = - 



{l-V 4 f) 2 dT 2 

(l + « 4 f) z 2 



/-, , 4 e .r 2 dy 2 + dx\ dz 2 
0- + V — ) 5 i + -¥32) 



CZ,T 



The horizon determined from g T r(v) — is at Vh = 
(^j-) 1 / 4 , thus it is moving away from z = (the AdS 
boundary) as needed. The 4-th power of v is related to 
the fact that its expansion near z — to the 4-th or- 
der is responsible for the stress tensor as observed on the 
boundary, which was tuned to correspond to the Bjorkcn 
boost invariant solution of ideal hydrodynamics [2 11 ] : the 
starting point for JP. 

This metric provides an asymptotic (large r) solution 
to the Einstein eqns 



(33) 



After this metric is substituted to the l.h.s. one finds that 
all terms of the "natural" order of magnitude 0(t~ 2 / 3 ) 
cancel out, with only the higher order terms remaining. 
More specifically, we found that only the terms ~ 
1/t 2 are present, with rather compact expressions such 
as 



r 2 T T7 



4 c 



(3 + v y 



(34) 



with c = (— ) t 

V e ' 



CT : y — ► y , x j_ — > cx± 
The resultant metric is 



ds 2 = — 



dr 2 



1 



r 4/3 



r 2 dy 2 



dx 2 ± 



dz 2 
z z 



The massless particle moves according to ds 2 
in JP metric is 



dz 



1 - 



0, which 



(40) 



1 



We have assumed that the particle always starts from 
outside the horizon: z < ts This EOM is solved numer- 
ically for different initial conditions. (From here on, we 
always use r = 10 as initial time for numerical solution, 
since the metric (|39|) is valid asymptoticly r » 1) 

To obtain the analytical form of the asymptotic behav- 
ior, we define: 



z 

4" 

T 3 



(41) 



r 2 T 



Av 2 



{Z + v){v-Zf 



(35) 



(36) 



Please note that those terms are not only subleading at 
large r but also are much simpler than all the terms which 



and the EOM becomes 

1 - It 1 MT 1 / 3 1 M 1 / 4 



vT+i 4 u 3/ 4 3 r 2 / 3 



(42) 



Note vl — ► 1 as t — ► 00. Assuming the second term 
dominates the first term on the RHS, we obtain the 
asymptotic form u 

assumption. In terms of z and r, we have: 



1 — ^ t 2 / 3 , which confirms our 
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FIG. 7: trajectories of massless particles, with initial z co- 
ordinates: z(10)=0.1(solid red) z(10)=0.5(dash-dotted blue) 
z(10)=0.9(dashed green) The horizon is also plotted (dotted 
black) for comparison. The trajectories of the massless par- 
ticles approach each other asymptotically, but does not seem 
to approach the moving horizon. 



Z = T 1 / 3 



1 - 



6V2 



.-2/3 



(43) 



For massive particle, the action is given by S = m J ds. 
Similarly we focus on the case that particle moves in a 
trajectory with constant y and x±: EOM follows from 
variation on action. Let z = r 1 / 3 /, then the function / 
needs to satisfy the following eqn: 



-27r 2 / 16 / 2 - 6r/ 17 / + 18r 2 / 8 / 2 - 108r 2 / 12 / 2 - 6/ 14 
-4/ 10 + 54r/ 5 / - 54r/ 13 / + 12r// + 108r 2 / 4 / 2 



-6r/ 9 / + 6f b - 3/ 18 + 9r*f"f + 9r z ff - 9r 
-18r 2 / 9 / - 126r 4 / 3 / 12 + 9r 4 / 3 / 20 + 27r 4 / 3 / 16 
-27r 4 / 3 / 4 + 126r 4 / 3 / 8 - f 2 + 9r 2 / 2 = (44) 

It is again solved numerically, with initial conditions 
rV3 a „ A ^ m^=. Note that 



r 4/3 



satisfying zq < t and z(tq) < 



l-i "n* 

1+ Tf) 4/3 



free falling massive object will move with speed of light 
asymptotically. We expect (|4"3"]) to be the asymptotic so- 
lution. By plugging (gSJ) in flS]), we get the RHS: |r" 4/3 , 
which tends to zero as r grows 

Furthermore, we compare numerical solution with the 
asymptotic solution in Fig[HJ The two solutions agree 
well at large r. This confirms (|43[) is the correct asymp- 
totic solution. 



FIG. 8: trajectory of massive particles starting with / = 0.8 
and / = 0. (solid red) at t — 10. The trajectory is indis- 
tinguishable from the asymptotic solution (dashed black) at 
t ~ 15 



To study the falling string, we first parameterize the 
string by z = z(t, y). Instead of solving it this form. We 
recall our experience with non-scaling solution in AdS 
space. At large enough r, the edge of the string will be 
less important, with most part of the string falling freely. 
Therefore we ignore the y dependence of z: z — z(t) 

Defining / = , The EOM follows straightforwardly 
from the Nambu-Goto action with the metric (f39"|) . It is 
a quite lengthy expression, which we choose not to show 
here. 

We expect the same asymptotic solution (|4"3"]) . By plug- 
ging gSJ) in the EOM, we get the RHS: -^t- 2 / 3 , 
which tends to zero as r grows. Fig[9] compares numeri- 
cal solution with the asymptotic solution, which confirms 
it is the correct asymptotic solution. 

Now we proceed to our final membrane falling 

in JP metric. Let zlr, y, x 2 , x 3 ) = r 1 / 3 /(r): the EOM is 
again quite lengthy and not shown here. 

We have solved it with a number of initial conditions 
and found that all extra terms are subleading near hori- 
zon, so this EOM gives the same asymptotic solution as 
the other cases, namely / = 1 — )t~ 2 / 3 . 

The numerical solutions are displayed in Fig llOl which 
confirm the asymptotic solution. 

We found that in all cases studied - massless and mas- 
sive particles, string and membranes - their late time 
behavior can be approximated by the same asymptotic 
solution 



(z- Zh ( T ))~[- 



6V2 



-1/3 



(45) 
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10.0 12.5 15.0 17.5 20.0 

time=tau 



FIG. 9: trajectory of string with initial condition / = 0.8 and 
/ = (solid red) at r = 10. The trajectory is indistinguishable 
from the asymptotic solution (dashed black) at r ~ 15 




10.0 12.5 15.0 17.5 20.0 

time=tau 



FIG. 10: trajectory of membrane with initial condition / = 
0.8 and / = 0(solid red) at r = 10. The trajectory is indis- 
tinguishable from the asymptotic solution (dashed black) at 
t ~ 15 

V. SUMMARY 

This is the first paper of the series, devoted to quan- 
titative formulation of the "gravity dual" to high energy 
collisions of macroscopically large bodies (heavy ions). In 
it we have formulated the setting in which the problem 
is simplified sufficiently to be solvable. 



Its central idea is that various "debris" from a col- 
lisions, in form of massless and massive particles or 
"stretching" open strings, all fall toward the AdS cen- 
ter. Although qualitatively such falling may look quite 
similar, the equations of motion and solutions are differ- 
ent for different objects. The main result of this work is 
a systematic demonstration of this statement in detail, 
both for initial time (when the underlying metric is sup- 
posed to be close to AdS) and at the late times (when 
the metric is close to JP solution) . As we will see in sub- 
sequent papers later, small differences in "falling" leads 
to quite different "holograms" in form of stress tensor at 
the boundary. 

One possible solution can be to unify all such "debris" 
as a single massive "membrane" , falling under its own 
weight. As shown first by Israel [33[ long ago, in such 
case one can greatly simplify the gravitational aspect of 
the problem, using two different solutions of the source- 
less Einstein equations inside both space-time domains, 
appropriately matched at the hypersurface made by the 
world-volume of the membrane. Two solutions are sub- 
ject to "junction conditions" providing new EOM for the 
membrane itself. We will discuss those issues elsewhere. 

Let us now point out few more specific results of this 
work. In the study of longitudinally stretched strings we 
have found that "scaling" solutions used previously for 
determination of "kink" 's anomalous dimensions are not 
at all adequate in Minkowski time. We found that while 
for wall rapidity Y > Y max sa 1/2 these solutions are 
absent, and there are two of them for smaller Y. We 
further studied stability of the solutions and have proven 
that at least for Y > Y c ~ 1/4 they indeed are unstable. 

Our main finding for generic non-scaling solutions 
(which come from numerical solutions of PDEs) is that 
while at small velocity of stretching there is the so called 
scaling solution, generically at high stretching one gets 
instead asymptotic approach to a "rectangular" solution, 
consisting basically of two near- vertical strings and freely 
falling horizontal part. 

Another result which was not expected is that although 
all types of objects - massless and massive particles as 
well as open strings and membranes - approach the JP 
horizon in the same universal way. Unlike in the textbook 
case of the Schwartzschild metric, this approach does not 
happen exponentially but only as a power r~ 2 / 3 of time. 
Note that this power is the same as appears in subleading 
terms, ignored by JP at late time. It remains a challenge 
to find an appropriate vacuum solution to Einstein equa- 
tion complementing the late-time JP metric. 
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Except that in Israel's problem of non-stretching black 
hole the horizon is stationary, while in our case it is mov- 
ing. 

As usual for a gravitational collapse, in a co-moving 
frame the horizon is not important and is crossed, which 
is not important for us to follow in this work. 
The reader may ask why we don't refer to conserved en- 
ergy, which will make this much simpler: the reason is 
the next section would not have this avenue open for us. 
The situation in which two charges move in the same 
direction is just a Lorentz boosted static solution: in this 
case a square root of v in the Lorentz factor is of course 
obvious. 

Similar problems have been encountered by previous 
studies of jet quenching, and another way to deal with 
them, proposed in Hertzog et al [Tlj ] . takes advantage of 
the re-parametrization invariance to fine tune the perfor- 
mance of PDE solver. 
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